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Abstract 

We present analytically the exact energy bound-states solutions of the 
Schrodinger equation in D-dimensions for a pseudoharmonic potential plus 
ring-shaped potential of the form V(r, 9) = D e (j^ — ^\ + | by means 
of the conventional Nikiforov-Uvarov method. We also give a clear recipe of 
how to obtain an explicit solution to the radial and angular parts of the wave 
functions in terms of orthogonal polynomials. The total energy of the system 
is different from the pseudoharmonic potential because of the contribution of 
the angular part. The general results obtained in this work can be reduced 
to the standard forms given in the literature. 
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I. INTRODUCTION 



The solution of the fundamental dynamical equations is an interesting phenomenon in many 
fields of physics and chemistry. To obtain the exact Estate solutions of the Schrodinger 
equation (SE) are possible only for a few potentials and hence approximation methods are 
used to obtain their solutions [1]. According to the Schrodinger formulation of quantum 
mechanics, a total wave function provides implicitly all relevant information about the be- 
haviour of a physical system. Hence if it is exactly solvable for a given potential, the wave 
function can describe such a system completely. Until now, many efforts have been made to 
solve the stationary SE with anharmonic potentials in two dimensions (2D) and three di- 
mensions (3-D) [2-6] with many applications to molecular and chemical physics. The study 
of the SE with these potentials provides us with insight into the physical problem under 
consideration. However, the study of SE with some of these potentials in the arbitrary di- 
mensions D is presented in (cf. Ref. [7] and the references therein). Furthermore, the study 
of the bound state processes is also fundamental to understanding molecular spectrum of a 
diatomic molecule in quantum mechanics [8]. 

The Harmonic oscillator [9,10] and H-atom (Coulombic) [9-11] problems have been thor- 
oughly studied in D-dimensional space quantum mechanics for any angular momentum £. 
These two problems are related together and hence the resulting second-order differential 
equation has the normalized orthogonal polynomial function solution. 

On the other hand, the pseudoharmonic potential may be used for the energy spectrum 
of linear and non-linear systems [12,13]. It is generally used for discussion of molecular 
vibrations. Additionally, this potential possesses advantages over the harmonic and leads 
to equally spaced energy levels. One of the advantages of the pseudoharmonic potential 
over the harmonic oscillator is that it can be treated exactly in three as well as in one 
dimension. Besides, the Pseudoharmonic and Mie-type potentials [12-14] are two exactly 
solvable potentials other than the Coulombic and anharmonic oscillator [9-11]. 

An exactly complete bound-state solutions of the 3D SE with pseudoharmonic potential, 
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anharmonic oscillator-like potential with inclusion of a centrifugal potential barrier does not 
complicate the solutions, which are available in closed form. Recently, the solution was also 
carried out by using orthogonal polynomial solution method and also by performing a proper 
transformation procedures [15]. Besides, the analytical solutions of the D-dimensional radial 
SE with some diatomic molecular potentials like pseudoharmonic [15] and modified Morse 
or Kratzer-Fues [16] potential are also solved by selecting a suitable ansatz to the wave 
function [17]. 

Chen and Dong [18] found a new ring-shaped (non-central) potential and obtained the 
exact solution of the SE for the Coulomb potential plus this new ring-shaped potential which 
has possible applications to ring-shaped organic molecules like cyclic polyenes and benzene. 
The complete exact energy bound-state solution and the corresponding wave functions of a 
class of non-central potentials [19] have been solved recently by means of the conventional 
Nikiforov-Uvarov (NU) method [19-25]. 

Recently, Cheng and Dai [26], proposed a new potential consisting from the modified 
Kratzer's potential [27] plus the new proposed ring-shaped potential in [18]. They have 
presented the energy eigenvalues for this proposed exactly-solvable non-central potential 
in 3-D-SE throughout the NU method. The two quantum systems solved by Refs [18,26] 
are closely relevant to each other as they deal with a Coulombic field interaction except 
for inclusion of a centrifugal potential barrier acts as a repulsive core which is for any 
arbitrary angular momentum t prevents collapse of the system in any dimensional space due 
to this additional perturbation to the original angular momentum barrier. In a very recent 
works [28,29], we have given a clear recipe of how to obtain analytically the exact energy 
eigenvalues and the corresponding normalized wave functions of the Schrodinger equation in 
.D-dimensions with the proposed modified Kratzer plus ring-shaped potential [28] and the 
modified Coulomb plus ring-shaped potential [29] by means of the conventional NU method 
[19-29]. 

The purpose of the paper is to solve the SE for the pseudoharmonic potential plus this 
new ring-shaped potential offered in [18]. This new proposed potential falls among in the 
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class of non-central potential which have already been solved in our previous work [19]. In 
spherical coordinates, we have given this appellation to the non-central potential 



where D e = kr^/8 is the dissociation energy between two atoms in a solid, r e is the equilib- 
rium intermolecular seperation and (3 is positive real constant. The potential (1) reduces to 
the pseudoharmonic potential in the limiting case of (3 — [13,15]. The NU method [19-29] 
has been used to solve the SE for this new potential (1). 

This work is organized as follows: in section II, we shall briefly introduce the basic 
concepts of the NU method. Section III is mainly devoted to the exact solution of the 
Schrodinger equation in D-dimensions for this quantum system by means of the NU method. 
Finally, the relevant results are discussed in section IV. 



The NU method is based on reducing the second-order differential equation to a gener- 
alized equation of hypergeometric type [20]. In this sense, the Schrodinger equation, after 
employing an appropriate coordinate transformation s = s(r), transforms to the following 
form: 



where a(s) and a(s) are polynomials, at most of second- degree, and r(s) is a first-degree 
polynomial. Using a wave function, ip n (s), of the simple ansatz: 




(1) 



II. BASIC CONCEPTS OF THE METHOD 




(2) 



VVi(s) = <f> n (s)y n (s), 



(3) 



reduces (2) into an equation of a hypergeometric type 



a(s)y'^s) + r(s)y' n (s) + Xy n (s) = 0, 



(4) 



where 
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r(s) = r(s) + 27r(s), r'(s) < 0, (6) 

and A is a parameter defined as 

A = \ n = -nr'(s) - V (s), n = 0, 1, 2, .... (7) 

The polynomial r(s) with the parameter s and prime factors show the differentials at first 
degree be negative. It is worthwhile to note that A or A n are obtained from a particular 
solution of the form y(s) = y n (s) which is a polynomial of degree n. Further, the other 
part y n (s) of the wave function (3) is the hypergeometric-type function whose polynomial 
solutions are given by Rodrigues relation 

R d n 

where B n is the normalization constant and the weight function p{s) must satisfy the con- 
dition [20] 

d t(s) 

—w(s) = —-w{s), w{s) = a(s)p(s). (9) 
ds ays) 

The function n and the parameter A are defined as 



, <T ' {s) - f{s) ) 2 -a( S ) + ka(s), (10) 



A = A; + 7r'(s). (11) 

In principle, since n(s) has to be a polynomial of degree at most one, the expression under 
the square root sign in (10) can be arranged to be the square of a polynomial of first degree 
[20]. This is possible only if its discriminant is zero. In this case, an equation for k is 
obtained. After solving this equation, the obtained values of k are substituted in (10). In 
addition, by comparing equations (7) and (11), we obtain the energy eigenvalues. 



III. EXACT SOLUTIONS OF THE QUANTUM SYSTEM WITH THE NU 

METHOD 



A. Seperating variables of the Schrodinger equation 

The proposed potential (1) can be simply rewritten in the form of isotropic harmonic 
oscillator plus inverse quadratic plus ring-shaped potential as 

V(r,9) = ar 2 + ^ 2 +(3^^- + c, (3 > (12) 

where a = D e r~ 2 , b = D e r 2 and c = —2D e . 

The potential in (12) solved for the limiting case of (3 = by using the orthogonal 
polynomial solution method [15] and by the Rydberg-Klein-Rees (RKR) procedures [13]. In 
fact the energy spectrum for the potential in (12) can be obtained directly by considering it 
as one case of the general non-central seperable potentials discussed previously in [19]. 

Our aim is to derive analytically the exact energy spectrum for a moving particle in the 
presence of a potential (12) in a very simple way. The D-dimensional space SE, in spherical 
polar coordinates, for the potential (12) is [28-30] 

fi 2 i 1 d f n ^ d\ 1 



— <— l/" 1 - I +- 

2/i [ r D 1 dr \ dr J r 2 



1 d 2 2/z/3 cos 2 9 



sine 06 y U d6) 
+ p (e - ar 2 - A - c) | V(r, 9, tp) = 0, (13) 



sin 2 9 dip 2 h sin 2 9 

where \i = ™^™£ 2 being the reduced mass of the two particles and ip(r, 9, (p) being the total 
wave function separated as follows 

^ n£m (r, 9, <p) = R(r)Y e m (9, <p), R(r) = r^ D ^/ 2 g{r), Y e m (9, <p) = H{9)${ip). (14) 

On substituting equation (14) into (13) leads to a set of second-order differential equations: 

^ + AW = «, (15) 
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1 d ( . d N 
s^9~dB V m ~dl 



m 



_ 2n(3 cos 2 9 
sin 2 9 ~~ Ti 2 sin 2 9 



+ £(£ + D - 2) 



# (0) = o, 



(16) 



1 d d\ £(£ + D — 2) 



r D 1 ^ r 



dr 



E — ar — c 



i?(r) = 0. (17) 



The solution in (15) is periodic and must satisfy the period boundary condition &((p + 2ir) 
$(y?) from which we obtain [9,10] 

1 



exp(±irrup)i m — 0,1,2, 



(18) 



Further, equation (16) representing the angular wave equation takes the simple form 



d 2 H(9) cos 9 dH(9) 



d9 2 



sin 9 d9 



+ D-2) 



m 2 + (2fip/h 2 ) cos 2 9 
sh?9 



H{9) = 0, 



(19) 



which will be solved in the following subsection. 



B. The solutions of the angular equation 

In order to apply NU method, we introduce a new variable s = cos 9, Eq. (19) is then 
rearranged as the universal associated-Legendre differential equation [26,28,29] 

d 2 H(s) 2s dH(s) u'(l-s 2 )-m' 2 



ds 2 



1 — s 2 ds 



+ 



sin 2 # 



-H{9) = 0, 



(20) 



where 



v 1 = £'(£' + D - 2) = 1(1 + D - 2) + 2n/3/tT and m = J m 2 + 2^ /H 



(21) 



The solution of this equation has already been solved by the NU method in [26,28,29]. 
However, the aim in this subsection is to solve with different parameters resulting from the 
.D-space-dimensions of Schrodinger equation. Upon letting D = 3, we can readily obtain the 
standard case given in [26]. By comparing Eqs. (20) and (2), the corresponding polynomials 
are obtained. In terms of the variable s, we have 



f(s) = -2s, a(s) = 1 - s 2 , a(s) = -us 2 + u 



m 



J2 



(22) 



Inserting the above expressions into equation (10) and taking a'(s) = —2s, one obtains the 
following function: 



7r(s) = ±yj(l/ - k)s 2 + k-u' + 



m 
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(23) 



Following the method, the polynomial n(s) is found to have the following four possible 
values: 



tt(s) 



(24) 



m's for k\ = v' — m' 2 , 
—m's for hi — v' — m' 2 , 
m' for ki — v'i 

—m! for k 2 = v' . 

Imposing the condition t'(s) < 0, for equation (6), one selects the following physically valid 
solutions: 



k\ — v' — m! and tt(s) = —m's, 



(25) 



which yields form equation (6) 



t(s) = -2(1 +m')s. 



(26) 



Using Eqs (7) and (11), the following expressions for A are obtained, respectively, 



A = A n = 2n(l + m) + n{n — 1), 



(27) 



A = v' - m'(l + m'). 



(28) 



We compare Eqs (27) and (28) and from the definition v' = £'(£' + D — 2), the new angular 
momentum £' values are obtained as 



e = _ {D 2) + _ — - — - + m/ + - 
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which can be easily reduced to the simple form 

£' = n + m', 
8 



(29) 



(30) 



in 3D [26]. Using Eqs (3)- (5) and (8)- (9), the wave function can be written as, 
H m i (6) 



(2£ f +!)(£'- m')\ 



\ ' 2(£' + m')\ 



sm m '(9)p( m '' m '\cos9), 



(31) 



where 



n 



(32) 



2 2 

with £' is given explicitly in (29). 

Hence, we are left to solve Eq. (17). It transforms, after lengthy but straightforward 
algebra, into the following simple form [30]: 



d 2 g{r) 
dr 2 



+ 



2fi 



-V(E-c) — 
h 2 V J h 2 



2fxa 2 v + (2fib/h 2 



g{r) = 0, 



(33) 



where 



v = -(M - 1){M - 3), M = D + 21. 



(34) 



Obviously, the two particles in Eq. (33) interacting via anharmonic oscillator potential plus 
inclusion of a centrifugal potential acts as a repulsive core which for any arbitrary t prevents 
collapse of the system in any space dimension due to this additional centrifugal potential 
barrier. Hence, Eq. (33) has to be solved by using the NU method in the next subsection. 



C. The solutions of the radial equation 

The aim of this subsection is to solve the problem with a different radial separation 
function g(r) in any arbitrary dimensions. We now study the bound-states (real) solution 
E > c of Eq. (33). By employing the suitable transformation, s = r 2 , and letting 

2/x6 



2a. 2ua 
s=J ¥ (E-c),a = i;r , 1 = u +1 _ T , 



one can transforms Eq. (33) into the following form 
d 2 g(s) , 1 dg(s) , 1 



(35) 



ds 2 ^ 2s ds 



(2sY 



i i 
e s — as — 7 



g(s) = 0. 



(36) 
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To apply the conventional NU-method, Eq. (36) is compared with (2), the corresponding 
polynomials are obtained 



t(s) = 1, a(s) = 2s, a(s) = e 2 s — as 2 — 7. 



(37) 



The polynomial ir(s) in Eq. (10) can be found by substituting Eq. (37) and taking a'(s) = 2. 
Hence, the polynomial ir(s) is 



tt(s) = i ± \^Aas 2 + A(2k - e 2 )s + 4 7 + 1. 



(38) 



According to this conventional method, the form of the expression under the square root in 
Eq. (38) is to be set equal to zero and solved for the two roots of k which can be readily 
obtained as 



(39) 



In view of that, we arrive at the following four possible functions of n(s) 



tt(s) 



\ + \ [2y/as + V47 + l] for fcx = f + \^Ja(k) + 1), 
\ - \ [2^s + v^TT 



for fci = e- + lJ a (4rf + l), 



2 

2 ' 2 ' 

for k 2 — y ~ iy/<*( 4r f + 1 ), 
2 



(40) 



for fc 2 = £ 2 



lv / «( 4 ^ + 1 )- 



The correct value of n(s) is chosen such that the function t(s) given by Eq. (6) will have 
negative derivative [20]. So we can select the physically valid solutions to be 



2~ 2 v ^ 7 ^ ^ = ^ as ' 



(41) 



which give 



t(s) = 2 + ^47 + 1 - 2v^s, r'(s) = -2y/a < 0. 



(42) 



Using Eqs (7) and (11), the following expressions for A are obtained, respectively, 



A = \ n = 2N^a, N = 0,1,2,..., 



(43) 
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e 2 1 / 



(44) 



Hence, the energy eigenvalues are: 



E N = c + , ^ Un + 2 + y/(M- l)(M-3) + 8/i6/^ 2 + 1 
^ 2/i V 

Substituting 

(M - 1)(M - 3) = Av = (D - 2f + 4£'(£' + D - 2) - 8/i/3/ft 2 - 1, 
with defined in (29), Eq. (45) transforms into the form 

E N = c + 



(45) 



(46) 



\ 2/x 



4iV + 2 + ^(,0 - 2) 2 + 4f (£' + D-2) + 8fi(b - f3)/h 2 , 



(47) 



where N,£' = 0, 1,2, ... Equation (47) contains the contributions coming from the angular- 
dependent part of the SE for the pseudoharmonic potential plus ring-shaped potential as 
well. 

(i) For 3-D-pseudoharmonic potential plus ring-shaped potential case, we use transfor- 
mation of parameters a = D e r~ 2 , b = D e r 2 and c = —2D e in Eq. (47): 



E N = -2D e + 



2h 2 D e 



2N+1 + ^(n + m' + 1/2) 2 + 2fi(D e r 2 e - [3)/h 2 ^J , (48) 



where N,n,m' = 0, 1, 2, .. and m! is defined in Eq. (21). 

(ii) For 3-D-pseudoharmonic potential case, we use the transformation of the parameters 
a = D e r~ 2 , b = D e r 2 , c = -2D e in Eq. (47): 



E N = -2D e + 



^ (2N + 1 + ^(e + 1) 2 + 2^D e ri/n 2 j , N,e = 0,1,2, 



(49) 



where £ = n + m and it is found to be consistent with Ref. [15]. 

In what follows, we attempt to find the radial wavefunctions for this potential. Using 
t(s), 7r(s) and a(s) in Eqs (5) and (9), we find the first part of the wave function 



4>(s) = exp 



a A s (nVJr+i)/4 



(50) 
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and the weight function 



p(s) = exp (-y/Eis) sv^ 1 / 2 



(51) 



which is useful for finding the second part of the wave function. Besides, we substitute Eq. 
(51) into the Rodrigues relation in Eq. (8) and obtain one of the wave functions in the form 



y n (s) = 5 n exp (yfa) s-^^ (^ + ^ TT/2 exp (-yfa)) , 



(52) 



where B n is a normalization constant. Hence, the wave function g{s) can be written in the 
form of the generalized Laguerre polynomials as 



g(r) 



where 



2 



\J(D - 2) 2 + At' {£' + D-2) + 8fi{b - (3)/h 2 - 1 



(53) 



(54) 



Finally, the radial wave functions of the Schrodinger equation are obtained 



R{r) = C N , L r L -^/ 2 e^ 



a 



(v^r 2 ), 



(55) 



where Cjv,l is the normalization constant to be determined below. Using the normalization 

oo 

condition, / it^r^r 0-1 ^ = 1, and the orthogonality relation of the generalized Laguerre 



polynomials, / z v+1 e z [L^(z)] dz = 
o 



(2n+ri+l)(n+ri)\ 



, we have 



C 



N,L 



2 (v2jm:/hr e ) L+3/2 N\ 

\ T(L + N + 3/2) 



(56) 



Therefore, we may express the normalized total wave functions as 



\ 



( x /2jlD' e /hr 



L+3/2 



(2f + !)(£' -m')\N\ 



+ m')\T(N + L + 3/2) 



.L-(D-3)/2 



x exp 



l lP < r 2 | L (L+L 2) 



\ 2ft 2 r 2 



^^r 2 ] sin m ' (9)Pi m '> m "> (cos 0) exp(±im V ) . (57) 



ft 2 r 2 
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On the other hand, the wave functions in (57) can be reduced to their standard forms in 
[15] as one sets (3 = 0. Therefore, we finally have 



ij){r,e,<p) = 



(y/2jnr e /hr e ) (2£+l)(£-m)\N\^ L _ (D 
\J 7i(£ + m)\(N + L + 1/2)! " 



(D-3)/2 



x exp 



l lP < r 2 | L (I- -2) 



\ 2h 2 r 2 e 



J N 



sin m (^)Pi m ' m) (cos^)exp(±imv9). (58) 



where 



L + \ = ^(£+l/2y + 2fiD e r 2 Jh 2 , i = n + m. 



(59) 



By choosing appropriate values of the parameters D e and r e , the desired exact energy spec- 
trum and the corresponding wave functions can be calculated for every case. 



IV. CONCLUSIONS 

We have easily obtained the exact bound state solutions of the D dimensional radial 
SE for a diatomic molecule with the pseudoharmonic potential plus ring-shaped potential 
by means of the conventional NU method for any N£'m'. The NU method can be applied 
systematically to both radial and angular parts of the wave function. Further, the wave 
functions are expressed in terms of special orthogonal functions such as Laguerre and Jacobi 
polynomials [31]. The presented procedure in this study is systematical and efficient in find- 
ing the exact energy spectra and corresponding wave functions of the Schrodinger equation 
for various diatomic molecules. This new proposed potential can be reduced to the standard 
pseudoharmonic potential in the 3D, which appears to describe the molecular vibrations 
quite well, by doing the following transformations: w! — > m, £' — > £ = n + m. We point out 
that this method is simple and promising in producing the exact bound state solution for 
energy states and wave functions for further anharmonic oscillator-type potential plus new 
ring-shaped potential. 
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